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I. — Introductory. 

The present paper originated in an attempt to discover the significance of 
^certain functions developed by Professor Karl Pearson in a memoir entitled 
"' A Mathematical Theory of Kandom Migration," and by him called 
W'fimctioQis, They belong to the category of normal functions, and are 
applied in the memoir named to obtain an expansion to represent a distribu- 
tion symmetrical about a point in a plane. The distribution is not fixed but 
depends on a parameter a, the function ftj2n being a function of ^ = rj g-, r 
being the distance from the centre of the distribution. The fundamental 
differential equation is {cP/dx^-i- (x-{-x~^)d/dx-{-2(7i-\"l)} ftJ2n (^) = 0. 

In the course of the present investigation it soon appeared that the same 
function led to solutions of the equation of conduction of heat in two 
dimensions for the case of symmetry round the origin, the time t taking the 
place of a^. In fact, it was found that if a solution of that equation is 
sought in the form /(t)(j) (r^/t), that solution is t~ (^+^) o)2n('^^/'^t), n being' 
arbitrary. The function (W2n is equal to g-^'Y^*^ multiplied by a polynomial, 
and is therefore especially adapted to the solution of the problem of the 
■cooling of an infinite sheet, the temperature at a great distance being always 
zero. 

The next step in the paper is to generalise the co-functions ; and all the 
solutions of the equation "^hi = du/dt are found which are of the form 
.f(t) ^ (r^/t) (y), where @ is a function of the angular co-ordinates of the point 



1908.] a Class of Normal Functions occurring in Statistics, 311 

alone. The result is that, as in other problems, @ is a spherical harmonic, 
or circular function, in the case of three and two dimensions respectively, 
that/(/()is of the form t~^, while (/> {r'^J4:t) satisfies the generalised equation, 

d^y/dx^ + {z-\-z~''^) dyjclz-\-{2n—m^z~^) y = 0, 
where z = Tt~'^ or xd^y / dx^ -]~ (x -{- 1 -^ m) dy j dx -{■ {n -{■ \m) y = 0, if « = r^^^- 
This equation is solved by the Laplace transformation, and the solution, 

(27rL)~'h~'^\e~^^ (i^ + 1)'^*+2"^~%2^^~M?^, taken round a loop contour from infinity 
round i^ =: 0, is denoted by con-i^m- The asymptotic expansion of con^m is 
{e-'^i-xy-i'^'/r 0? — im+ 1)} 

.. __ '?^^ — {^my , { ^^ — (i^^'O^ } {(^""l)^"? ( 2 ^^^ )^ } _ 
j_ . -J- , . _^ ^ ^ ^ 

which terminates if either n—^m- or n-\-\m is a positive integer. 

Thus the solutions of v% = 'buj'dt in the cases of one, two, and three 
dimensions respectively are 

r (-+%ia),+^, J {x\ r (-+%^-a),, .. ix) "^^^ 77^0, r (--^%i^a).,_ X, ,^ . ix) 7, (^, </>), 

sm 

^u being a spherical harmonic of order /v, which reduce for certain values 
of the constants to simple well-known solutions, viz. : — 

(i) To {-ivtyh-'^iioY n = ~-i), (ii) to r^e"^ (for n = m =0), 

(iii) to t'^eT"^ (for n = ^, ]c = 0). 

The one-dimensional solution obtained is not new, it practically coincides 
with the function of the parabolic cylinder, or the functions considered by 
Sturm* in connection with the conduction of heat, by Thiele and Charlier 
in reference to probability, and by Hermite.f The two-dimensional functions 
seem to be new, but they are shown to be the equivalent in polar 
co-ordinates of the functions considered by Hermite 

and the relations connecting the two types of function are obtained. 

It is shown in the paper that any function of x and 6 can be expanded in 
the form 

2S(X,,^ „^(t>.,^, ,j^ cos m6, where ( — ) ira,^^ m= \ /^^«~*^''a),,, _„i cos m6 dx d9, 

and a proof of the convergence and validity of the expansion is given in the 
case where /is independent of 0, with an indication of the extension to the 
general case. Hence a series is obtained giving a solution of the equation 

^ Liouville, I and II. 

t ' Comptes Eendus,' vol. 58. 
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y% = du/di, which reduces for t = to to Si function given all over the plane,, 
viz. : — 

XXcCn^ m (^o/O ^^^ ^n, m(^'^/44) COS mff, the coefficient a.y^^ 5^ being obtained as 
above. This series is shown, by direct expansion, to be equivalent to 
Laplace's definite integral solution of the cooling of an infinite plate, viz., 

\f(^, V) e-^(^-^)'+ (v--r,)'VHt-to) cl^drj/iTT (t--fo). 

Various other expansions are given analogous to JSTeumann's expansions of 
Jm (r^-hp^—2rp cos 0)^ in a series of Bessel functions. 

The last part of the paper considers the equivalent functions in Cartesian 
co-ordinates, in terms of which Hermite has expanded any function.. 
It is shown that an equivalent method of obtaining Hermite's expansion 
is that of equating the successive moments of the series to those of the given 
function, with special reference to obtaining the best fitting surface for 
a given frequency distribution in two dimensions. The series of approxima- 
tions so obtained is compared with that of Edgeworth,* which is ultimately 
the same, and is shown to converge rapidly under the same conditions. An 
example is then given of the contour lines of a surface obtained by taking, 
terms up to the fourth order. These contour lines show much similarity 
with those actually observed in statistical work, such as those given by 
Perozzo in his memoir.f The possibility of obtaining a fairly close approxi- 
mation to a given distribution by means of a knowledge of the product, 
moments up to the fourth order is thus established. 

II. — -The Derivation of the Generalised m-ftinction, 
1. The Fufhdamentcd Diffeo^ential Uqitation. — In order to obtain the types- 
of function suitable to the problems mentioned, the most general solution 
of the equation "^ho = dti/dt will be found in the form T.E.B, where. 
T is a function of t alone, E a function of r^/t, and @ of the angular 
co-ordinates. To begin with, the cases of one, two, and three dimensions are 
considered separately. 

(a) One Dimensio7i,—'Pxiiting z = r^/t in the equation, it becomes 

If to ^ Z,T, Z and T being functions of z and t respectively, we have 

so that both sides must be independent of t and z. 

^ ' Camb. Phil. Trans.,' August, 1904. 
t ^ Annali di Statistica,' 3a, 5, 1883. 
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Hence t~-{-nT =^ and 4tz~ +(z + 2)^-^nZ = 0, 

at dz^ az 

where oi is any constant. Thus T = t~'^. 

(b) Two Dimensions. — The equation in polar co-ordinates is 

3% 1 die 1 3% _ du 
3^2 r dr r^ d6^ dt 

lito = @v satisfies this equation, @ heing a function of 6 only, and v heing 
independent of 6, the functions B and v satisfy 

m being a constant. Thus © = A sin m^ + B cos m^^ and again, putting 
z == r^/t, and v = T .Z, the same procedure as above gives T = t"^ and 

(c) Three Dimensions. — Exactly similar procedure in this case leads to the 
conclusion that, if ?^ = T . Z . ©, @ being a function of 6 and 0, must be a 
spherical harmonic — say of order h, T must be t~^, and that then Z satisfies 
the equation 

(d) Three Dimensions — Cylindrical Co-ordinates (p, 0, ^). — In this system 

of co-ordinates the general solution is mO , Jc^Z, where Z satisfies the 
equation 

The equations (1) to (4) can be brought under the same form. Putting 

Z = z^X, 2;^^*X, #X, z^'^X in the four cases respectively, and also z = 4:X, 
the equations become, using dashes to denote differential coefiicients, 

xX'' + (^ + 3 ) X' + (n 4- J) X = 0, (5) 

xX'' + (x -{- 1 + m) X' -f- (7^ -f- ^m) X = 0, (6) 

xX''-{'{x + ^-i-k)X' + (n-}-ik)X = 0, (7) 

^X''H-(^+l + m)X'+(^ + P + im)X = 0. (8) 

Taking (6) as the standard form, the particular solution dealt with in this 
paper will be denoted by con~i,m' the corresponding solutions of (5), (7), and 
(8) are therefore 
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and the complete solutions of the original equations are 



oXXi 



(9) r (-+%*a),+., J {ml (10) r (»+%4- 21 '^^ «% ^(Pl 

(11) r <-+%i^a),^j,,^| {x) yu (0), (12) r C-+%i-«.,+,.,,, (x) ""^^l hz J?^ m^. 

The simplest solutions are given in the four cases as follows : — 
(i) n = — I, (ii) n = 0, m == 0, (iii) n z= ^^h = (), (iv) % = m = Ij = 0. 

2. .5%e SokifAon of the Funckimmital MqtmHon,'-~~'Puttmg X = e~"^^'^dv 

J c 

and forming the Laplace transformation of the equation, we find that V must 
satisfy the equation 

V (t?— 1) dY/dv = {(m— 1) -y— (^^4-|m— 1)} V, 

giving T = A-y^+i'^^-i (v— l)i«*"% 

while the contour c must be such that Ye'~^%V'—v^) has the same value at the 
two extremities. Assuming x to be real and positive, this will be attained by 
taking a contour consisting of a line following the real axis from i? = oo to 
v = a(a<l), then a circle about v = l of radius a, and a straight line 
returning along the real axis to infinity. A similar contour encircling the 
point instead of 1 will be denoted by 7. Putting ttz=zv—l and A =: ^iri, 

we have for mn^r^ the expression (^irije^^l 0""*%^^ + l)^'^i^iti''''~n~MtL 

J To 



The integral gives the asymptotic expansion for Wn,m 

which terminates if either n—^m or n+^m is a positive integer 

The most important case is that in which ^ — Im is a positive integer; in 
that case 

.^j 1 Q) 

(14) mnm = e ^r(^+tm + l)-< .-^r — . .x/ . — 7?— iT? oTTtt t tti 

^ ^ ' ^ ^ {T{m+l)(n—^m)l r(m + 2)11 (-/?.-. |m—l)! 

r('W.4--|^ + l)(^— |m) !. 

In particular, if m = this reduces to Pearson's function Q)2n* 

When 01— ^m is half an odd integer, m^^r/k becomes purely imaginary. This 

case occurs below, but no difficulty arises on that account. 

Especially important is the case in which both n^—^m and m + |m'"are 

integers, for then we have the relation, 

(15) {-oj^^'''o)n,jr(%+im + l) = (^x)-^^^'wn^^JT(n-^im + l), 
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3. Relations connecting Successive Functions, — 
{a)DiffeTence equations— 



tl 



y 

= 27n {(72-fim) «^,,H — 0^ — i^?0 ^n-l,m}' 
Thus X do)n-i^ ^jclx + (ll — |7?^) COn- 1, m = 0^ + 1^^^) (*^n. m- (1 ^) 

Applying this twice and simplifying by means of the fundamental 
differential equation, we find 

(7^-f -|m) 0)n^,n—(2n — X—l) «n-l,,;^+(7^ — |m— 1) C0n-2,m = 0. (17> 

(h) The Expression of (On,m <^s a Differential Goefficient. — Since the 
derivative of any solution of v^('^0 == 'ditjdt with respect to t is also- 
a solution, and since djdt {t^f(r'^]t) } = f^'^i^ O'VO? we deduce at once, 
after adjusting the numerical coefHcients, the following relations in the cases 
of one, two, and three dimensions respectively (^vicle (9), (10), (11) ) — 

a/a^ { r ^^^0)^-1, X (^)} = -7^r ^-i^W+i J (^), (18) 

a/a^ {r V«^o)^_i,^(^)} = ^{n-\m) r^- V^o).^,,^ {x), (19) 

a/a^ {r«^^^o)^_.,,+^(^)} = ~(7i-p) r^-i^^V-i/.+i (^). (20>= 

III. — :7%^ Linear Function w^+i,!. 

1. Eestricting the work in this section to the functions arising in the 
one-dimensional case, which are briefly called the linear w-functions, it wilL 
be convenient to write w^ for w^-i^i and (d\ for ft).^_i^_i,.and v in place of n. 

From (18) we then have 

t---ix^co, = (-d/dty {t'^x^coo (^)} r (-|)/r (z^+D 

= (-y-Hd/dty(t-k-^)/v (i;+i) = (-)^-^(a/ar)2-(^5-^6-^^^/^o. 

assuming that v is an integer. If v — ^ is an integer, we have, similarly^ 

t-'-wco, = (•^y-h(d/dty-h{t-^x-i(Di(x)}/v(v+^2) 

= ( -)--^ (3/37^)2^-1 {irtr6-^V4i} ^ (^y-i(d/dry' (t-h-'^'^^'). 
Thus in either case 

0,., = (-.)v-i^i- (a/ar)- (e-^-^/^o/^^r {j (z^+ 1)}. (21) 

Thus, as a function of r, wi.. is identified with the function considered bv 

''2 V 

Hermite.* 

■^ * Comptes Eendiis,' vol. 68, 
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2. An Uxjpansion in a Series of Linear co-Functions. — The expansion of 
"Wy is --~^ . . ,,. i 1— !i^ Z-f^iiZ ^Y ^^^"^ .),,, L which, if v is 

■an odd integer, becomes 

vl f 1 4^; (-.4^)1(^-1 )- 

.and if v is an even integer, 

_^ z^!(— ^)~^ f 1 4:X _(i^\ 

' 2^r{-H^+l)}l>l i(z.-2)!2!"^-- ' vl J' 

Hence, a)|^ = coefficient of ^^ in z/ !e^^'+^^^^^-^/(— a?)^2^r {i(i/+l)}> 
■and therefore 

(22) 
^nd 

(23) 

the former summation being for even values of v and the latter for odd, or 
together, 

00 

^-l^(g+c^P ^ \/-7rX(#ft>i./r(iz.+ l). (24) 

v = 

3. Application of the above Ux2Mnsion.— The foregoing properties will now 
be applied to obtain an expansion of Laplace's solution of the problem of the 
'icooling of an infinite bar. We have 



{27r (t-tojyh-^''-^^^^^'-''^ = (27r)-i 



.00 



^- w2 (t- to) - M {r-p)gq^ 



GO 



.00 



= (47r) ~i i e-^'i-^^^peP'^^'o t^iHo^'-^^^coi, {p^l^h) dujV (|z^ + 1), 

J — oo 

aQO 



— 00 



= (-2to)-ie^ &xi % {to/ty+ico'i^ (^) «i, (x), 
iusing (15) and putting ^ for p^jUi^. 



.00 



IsTow, as Laplace as shown, V = /(^) {4t7rt)~h--^''~^'>^l^^d^ is a function of x 

J —00 

satisfying dW/dx^ = 3V/3^ and reducing when ^ = to f(x). Putting (^—^0) 
in place of t, and expanding the subject of integration as above, we obtain 
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the following series to represent a solutioii[of the same equation reducing 

when t = ^0, to/(^) : 

ta,{hlt)^^''-^^)xim{x\ (25) 



CO 



where a, = — J^ f{p)e^_^ia>\ (|) d^. 

J —00 

In accordance with a remark made above, it may be noted that, when co^ is 
real, cd\ is imaginary, so that the expansion is wholly real. A different 
method of establishing the same expansion is given below. 

IV. — The General CD-Ftmction — Integral Properties, 
Let 



-—^TT lytj fi'^ 9JJ = ^tTT I e (0^,1^ ,jj \Xj G)^/^ _5^ yx) ax 

Jo 



dx 



•' Y w y 

The asymptotic expansion of the w-f unction shows that the integration in 
regard to x to the infinite upper limit is valid, and the order of integration 
may be inverted. Thus 

j du I dv ('^6 4.ly'^+l»^'^^l»^-«-l(t; + ly'^'-J«^'y-i"^-*^'-l (-^^4.^4. l)-i^ 

Jy Jy 

For every value of v the subject of integration is zero to the second order 
at least when %t is infinite if m is positive, and vice versa if m is negative. 
Taking m positive for the sake of argument, 

~47r2I.,.',. 



' '*'7r J-}^J n'^ 1% — I aw 



m 



du dv (U + ly+h^Hm''-'^-'^-^ (v + l)^'-|m^.-im-.i'»l (^^^ 4. ^ .^ ;l^^|-I^ 
-v~l J y 

= (-y^27rt f (i'y 'y^-'^'-i (1 + i;)^'-^-"i = (—)»^ (-.477-2), 

Jy 

if n = t?.', or = if ^^ ± n\ Thus 

^""^W^t, m (^) «Wb', -m {^) ^^ = 0, (26) 

Jo 




»co 



/•CO 

but e'^COn, m (^) W«, -m (^) C?^ = ( — 1)^. (27) 

Jo 
Again, consider the moment integral 

}C0 <»CO /« 

Jo Jy 

== r (r + 1) [ {u-\-Vf'-^"^'^'-''~H^'''-'''~^ du. 

J V 



YUJj. iiXXXl. A» 
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Thus, if n-\-\m—T is a positive integer, (l-j-i^^w+lw*-;--! ig a polynomial and 
the subject of integration has no residue at ti = 0. But if oi + ^m—r is zero or 
a negative integer, the residue is (— )''~^"'(?'+l— n — Jm)...(7'— m)/(7i— |m)!. 
Thus 

T(r-n-im + l), (28) 

according as r <, =,or > n-j-^m. In this r is not necessarily an integer, 
but n—^m and 7^ + ^m— r are assumed to be integers. 

If these results, (26) and (27), be applied in the one-dimensional case 
considered above, assuming that /(f) = 2(^^f*a)j,(|), we obtain at once 

a^=z^iii /(|)e^Ma)'^(f)^f. Consequently, if ^ z= r^/4:tQ and .v =: r^/4:t, 

Jo 

Xcc^{to/ty^'^'^x^Q)v{x) is a function satisfying d^u/di^^ = dic/di, and reducing 
when t = to to the given function /(^). Thus the expansion (25) is obtained. 

V. — The m-Functions in a Plane, 

Coming to the special values of n and m suitable to the development of 
functions in a plane, m is an integer and n—^m is an integer. Using 
relation (19), we have 

=: {-djdty {ri (rVO'^'^^o.m^"^^^} r (|m+ l)/r (Im^n^l) 
= { - V^Y{r^ (rytyi^^'coo, ,^e--'} r (|-m + l)/r ( J7?i + 71 4- 1). (29) 

In particular, if we start from the function independent of 0, i.e. put m = 0, 
O,,,o = r«+ia)« (ryU) = i-d/dtf {t~h-'^'t^') j n\, so that 

mn={-tY/nl(dydr^ + l/rd/drfe~^'^^K (30) 

This particular set of functions may conveniently be called the zonal 
«()-f unctions. 

If m be not zero, but n—^m a positive integer, so that a)^^»^'^ is polynomial, 

== (_a/a^)n-|»« ^^.m^-m-l-"rV4i^m0j. r(m+ l)/r (|-m + 7^+ 1) 

=: ( _ )n-hm (y2)^-im ^^„m^-m-l^-rV4f^tm0 j. p (^^^ ^ 1)/ p (1.^^^. + 9^+1). (31) 

Thus all the functions found are expressed in terms of the derivatives of 
the fundamental function with respect to r, 6, and i 

Since the derivatives of any solution of v^'^^ = dujdt with respect to 
Cartesian co-ordinates x and 2/ are also solutions, we may expect all solutions 
to be expressible in terms of the derivatives with respect to x and y of the 
fundamental solution. 
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The following equations may be shown to hold for the planar w-f unctions, 

and, therefore, 71— \m being an integer, 

or {d/dx-\- id/dyy-^'^'^' (d/dx-^-id/dijy-^^^^ (r i<?-^"/^0 ; ('^2) 

e.g., we have 

{d^ldx^-d^jdif) {rh-'''/^') = (4^J2)-i {^^2ny-TVu COS 29, 

2 d^ldxdy (r i^-'-'/^O = (4^^)"^ {r^h) e"'"'^' sin 26, 
and (dydaP + ^/dif) {rh-'^'l^') = ---{f)-^ (l^r^jU) e-'^'l^\ 

giving the values of d^/dx^, ^^l^y^, d^/dz^ in terms of Xlio and II12. Thus the 
connection is established between the planar w-functions and Hermite's 
functions o'^'^^/dx'^dy^ j^ <3 -«(•»'+ 2/^) |^ to which reference is made below. 

YI. — The Exidctnsion of any Function in Planar ay- Functions. 
1. Suppose that a function of r and 9 may be expanded in the form 

/ (r, 9) ISanm (^)*"* (»ynm (pc) COS m9 + 'tthn, mJ^^'^'cOn, m (^) siu m^, 



n m n m 



where x = r^j^tt, t being a parameter. Then, using the Fourier method for all 
values of m, we have 

f{r, 9) COSm9cl9 = TrS^nm^^^^W^i, m 

Jo % ' 

and f(r, 9) sin m9 cl9 = irt'bn,,npc^'^ci)n, m. 

Jo 

In order to find the coefficients a^^^^^^ K,m, we now use the integrals 
developed above, and obtain 



( — y^ TTCCn^ ,n= if. C^X'^^^COn, -m COS m9 dx d9, (33) 

Jo Jo 
with a similar expression for hn^m- Thus, in order to obtain the coefficient 
of (*>n,mi we make use of the allied function con^-^my which bears to it the 
relation (15). 

The values of m and n required to be taken do not appear, but it will be 
shown presently that m takes all integer values and n takes those values for 
which n—\m is zero or a positive integer. The expansion is therefore one 
in a series of polynomials. The corresponding expansion in Cartesian 
co-ordinates is given by Hermite.* 

* * Comptes Eendus,' vol. 68. 

Y 2 
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2. The Expansion of the Xajplace Definite Integral in a Plane, — An expan- 
sion analogous to (26) will first be obtained, x standing, as everywhere, for 
r^l4:t, let z stand for r<3''^/2 \/t and z for Te~''^l2\/t. Then, taking all values 
of n and m for which n—\ni and n + ^m are both positive integers, 

con, m = e-"^ {n + -|m) ! \ -— -— - - - — f — — - ... ■ 



{n + Jm)! (7^ — -|-m) ! 
= (7^ + |-m) ! a^-*"* . coefficient of 7i^~^''Vi;^+^«^ in e^^+(A-A)V2:-a:. 

or t t x^'^'con^ .Ji^- W^^+^^7(^ + |m) ! = 6^^+ ^^) (^- ^^), (34) 

w = m = — 2n 

7^ taking all half -integer values from to infinity. 
Again, taking the same range of values for m and n, 

Expanding the last expression by Taylor's theorem, we have 

xi'^'e'^^'^conm = ( — )^+^*'^ (3/3^)^-^"* {d/dzy-^"^'"' {e-''')/(^i — Jm) !. 
Thus, putting ^ = re'^ and ?; = re'^, 

xi^e'^'^cDnm = ( — )^+^«^4^2^^+i (3/3^)^^-^*'* (3/ 3^)^+*^'^ {t-h"'''!^^)! {n — \m) L 
This gives the proof of the relations given above (32). 

Putting ih = p + iq and —ih^p—tq in (35), that expression becomes 
^'^■¥q^+i{pu-\-qv)-x^ where (^6'^;) are Cartesian co-ordinates. 

Hence 

fee >>(X) 

dp dq e~{p''-^cP){t-k)-''P{u-u{}-tq{v-v,} 
— co«' — CO 

= (27r)~i i rd'p dq e-(;p'+^') i-i(2?u+<?t^) g^i2:S^Ci^^''V'^^'^a)^^(^^i) h''-i'''J{f+i'''/ {n -f |-m) r 



-co J —CO 



by (35), where % = (^^iH %^)/4?^o. Now put h' = h^/itlto) and //= k^y{t-~tQ) 
and V = jp-f-i^, 7/ = p — cq. Then the above becomes 

:^:g (_)n-jm (3/3.)^-i^^(3/3^')"'"^"'(27r)-i 

foo /.CO 
^-(^2+^2) t-<px-\-qy)^lp g(^ t^H-'^-H'''^X^'^^e''^'^Wnm (^l)/(^^ + Jm) j 
— CO J —CO 

The summation throughout is for all values of m and n such that n-k-^m 
and ^— Jm are positive integers. There will thus correspond to each term in 
the summation another with an equal but opposite value of m, and, using (17) 
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these terms will only differ by 6*»'(<^~^) being changed to ^~»**(<^-^). Eemeni- 
bering also that n—^m is an integer, ( — )*"2» == (^^y~m^ Hence 

CO 

+ 2tt (-)Vixi^xr>G>nm («) Ci>t.,^m (%) COS m ((j>-'d) t^H"^^-^. (36) 

This expansion is suggestive of Neumann's expansion in the theory of the 
Bessel function. Another theorem is given below, which, in some respects, is 
a nearer analogue. The development of the Laplace definite integral follows 
at once from (36) and is identical with that in the last paragraph (33). The 
method of this paragraph shows that the values of m and n are as there 
stated. 

3. The Expansion of m,n,n{{r^—2rpm%{<l>--d)-^p^)}j4:t in. a Series of 
(O'Fiinctions, — As in the last section, 



*CO /.QO 



J — 00 J — '00 



.00 /.CO 



(the summation being taken for the same values of m and n as above) 



= (2^)-i 



,<» 



,00 



— 00 



dp dq (3/3|)^~**^* (3/8?;)'^+*»*^-Ci''+?') t-'t(pu+qv) 

p2n^mc^l(^^ — |m) ! (n + 1 m) ! 
= tt (— )^+*^4-'^r^^%2»*ft)^^^^gim(.A-0y^y(^_i^^)^^ 

= tt (-y-^^n~^pyM)--^x-^^^mn,n,(x)QOBm((^'^d)/(n+^^^^ 

where the summation extends only to values of ?i such that n—^m is a 
positive integer. Thus wojo {(^'^—2rp cos (^—^) + p2)/4^} is expanded in a 
series of ca-functions in the variable {r^/4:t), and from this may be derived, by 
means of relations (32), the expansions of (On,^i {(r^ — 2r p co^ ((j^ -' 0) + p^)/4:t}, 
using differentiations with regard to t, x^ and 3/1. 

4. The Convergence and Validity of the Expansion of a Ftmction in a Series 
of w-functions. — The proof of the validity of expansions in a series of normal 
functions is always a difficult matter, and a general theorem to cover them 
has yet to be established. A discussion is here given of the expansion in a 
series of zonal ca-f unctions of a function symmetrical round a point, which 
seems capable of extension to the general ca-functions in a plane without 
much increase of difficulty, though with considerable addition to the analysis. 
The multiplication of different discussions for each particular set of normal 
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functions seems undesirable^ however, and the present investigation will 
therefore be limited to the single case of the zonal functions. A different 
method appears possible in the discussion of the validity of the processes in 
the development of the Laplace integral given above. 

In the following, the basis of the work is the definite integral expression 

r--a;(l+w) 

for the ce)-function, viz., (2Tn)~^ e{U'^l)Ho~^^~'^dii>, 

Jy 

The expansion to be justified is 

00 c'^ 

"^ {x) =: ^arCOr {00), wlicrc a^ = f(y)a^ri'^)^^d^J' 



JO 

Consider the sum of (7^+1) terms of the series 

00 n r<=^ 

J 

= — (47r2)"-i e^fiy) dy du dv g"~(«+i)y-(«+i) x 



r--U,~r— 1 



t {^l^-ly {v-VlY ir'-^H 



-00 t C 

= — (477-2)-i ^~^f{y)dy\ du\ dv{ii-\'V-\-l)''^ 

Jo y y 

A necessary condition is at once obvious for the validity of the processes 
involved. For, if y is very large, since the 'i^-contour includes points for which 
the real part of ti is negative, e'^^y is not everywhere finite on that contour. 
But iifiy) be such that, no matter how great y may he,f{y)<e^y, where Jc 
is some finite negative constant, we may take the it-contour such that the 
real part of u is always greater than k, and then the integrations required are 
justified. 

Let 

J^(x,y)= [ \ e-'^y-'"'^ (to+iy (v-tiy u~-'v-'^ (to + v + l)--^ dudv. 

Then, if y is finite, Jo is zero, since the it-contour includes no point at which 
(u + v^l) vanishes for any value of v on the contour 7. Also, by the 
condition assumed to be satisfied by f{y), the part of lo arising from values of 
y greater than a finite magnitude rj can be made as small as we please by 
taking rj sufficiently large. Thus To = 0. In treating I^+i it has to be 
remembered that the limit of the integral is sought when n tends to infinity, 
and that since the contour 7 encircles the point 0, and the two contours for 
u and V have been assumed such that u-i-v+l does not vanish for a pair of 
values within these contours, | (u+l)/u | and | ('^ + 1)/'^ | are both greater 
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than unity on the major portions of the contours, so that the subject of 
integration therefore becomes indefinitely great as n increases. 

/•oo 

Consider now the integral J^ (x, y) cly and divide into two portions for 

Jo 

which y is greater and less than x respectively, calling them Kj and K^* 
Let Ka be the same integral with the limits a and x (a<x). 
Then 



K 



dpe'^ { clu [dv 6-«-«y {u + Vf (»+ 1)'' tr^'V'' (m+ v+ 1) "^ 



a 



r r 
— cU dv 6"^ (!+«+*') (u + iy (v + iy ^r^-i^-*' (ti+ v+ i)-v 



y "y 



+ [ du j dv e'^'^^^-^''^"^'' (tc+iy {v + lYu-~''-^H-~'' (u + v + l)^^^ 

= — iK«-i-2Ka (say). 

Consider the latter part 2Ka; it is absolutely and uniformly convergent, 
and may therefore be written 

--d/dx { dtc.dv 6"^<i+^>~^^*(t6 + lyC'y + iy'~ho-''+h-%u + v + 1)~^ = —dLr/dx 

(say). Within the contours considered there is no pair of points for which 
(u-i-v-i- 1) = 0. If the t6-contour be extended so as to include all the points at 
which {^c-{^v + l) vanishes for points on the t-contour, we add to L,. the integral 

j I ^~^a+«)-«w(^^4.1)»•'^r*'-l(t;+l)*'~^^~''('^^4-'y+l)~^^^i^f^'y, th^ contour 7, 

J yJ y, 

including to •= — 1;— 1, but not tt = 0, and this addition 



= —27n 



e~(^~a)i^+^)(l-^v)-hlv = 0, 



y 



since the contour y does not include the point ( — 1), and since x-^a is 
positive. Thus the '2^-contour in L^ may be supposed extended so as to 
include the point at which (u-\-v-\-l) vanishes for all values of v or 7. 

Let the -y-contour be now specified more precisely as that part of a circle of 
large radius K, with centre at -y = — ^, that lies on the positive side of the 
straight line on which the real part of -y = --^, together with that part of this 
line which lies within that circle. This contour is considered in three parts, 
a, /3, 7, as follows : a, the straight portion ; ^, that part of the curved portion 
adjacent to cc and subtending a small angle S at the centre of the circle; 
7, the remainder. 

For values of v on a, the tt-contour, extended to include — (l + 'y), is taken 
to be a contour just surrounding the complete -y-contour. For values of v 
on ^ and 7 this contour is extended in the negative direction so as always 
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to be a short distance on the negative side oi u = — -y— 1. The ^-contour 
is considered in three portions similar to a, /3, 7, viz. : ai, the straight portion ; 
y8i, the arcs contiguous to ai and subtending an angle S at the middle point 
of cci ; 7i, the remainder. Thus the whole of the double integral is divided 
into nine portions obtained by combining one of the partial contours a, /3, 7 
with one of the partial contours ai, /Si, 71. 

Now e~'^^^(l-\-v~'^y remains always finite as r and v become large, except 
for V •= —1, a value which lies outside the region under consideration. Thus 
|6-^^(l+^"iyi < AI^-^^+^/^l and l^-^^ (l + ^~^y| < B|6-«^+^/**l. Thus, if 
V = E6^^|e-^^(l + ^"^y| < Ae(-^»+*'/^)«««'^, and, it u = EV*^, \e-''''{l + ^^^y\ 
<Bg(-«E'+r/E')cos,^. Hence, if 

and therefore, on the arc 7, cos ^ being > S, 

|6-««(1 + ^""^)^| < B6-(«-^)^'^ < B6-(«-^)p«~^ (where E' = pB-^\ 

and this quantity tends to zero as S diminishes, provided p remains finite. 

On the arc /3i, the expression certainly remains finite. 

On the straight part ai, [ (w + 1)/'^ | takes its greatest value at the point 
where it adjoins /3i, ix, {(u-\-l)luY is finite on ai, but along a length (^E')* of 
i:t measured from the real axis it becomes comparable with (1 — 1/Z;E')^ ^.6. 
with {l — ljk'R'y^'^ i.e. with 6-^^', where E' is large and c finite. Similar 
considerations apply to the various portions of the '^-contour, and putting 
them together we deduce that each of the portions of the double integral 
tends to zero as r tends to infinity. Thus Ur^^Jur = for all values of x, 
and therefore s^a = —dLr/Bx also tends to zero as r becomes infinite, for any 
value of a less than x, and in particular if a = 0. 

Thus 

Ka — II ^-^(i+^>-^y (^6+iy (v+iyir''v-''(tc + v + l)-^duchdy 

•' aJ y*' "yi 

y y 

and is independent of the lower limit a. 

The u and 'y-contours being the same, this last integral will be unaltered if 
u and V be interchano^ed. Hence 

vr — _i 



1 {U + dRldx], 



r r 
1 

2 
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where 
H 



JyJyj 

where e tends to zero as r becomes large ; since, by reasoning similar to that 
used above, the double integral H with the 'y-contour extended to include the 
point —{1+to) vanishes as r tends to infinity. . 
Thus 

L^^ = ooH= --27rt —u~'^(u-i-l)~^du= --47r^ and L^y = oo9H/3^= 0. 

Jy 

Finally, therefore, Ltr = ^Ka = — 27r^, whatever the value of a. 
Apply now the same treatment to the integral K&, viz. : — 

me-xii+v)-uj/ (^ ^ ly (^ ^ ly ^;,-r^-r (u-i-V'{-l)du dv dy 



yy 



— j j ^-^ (i+t^)-6w (iiJ\^ ly {v-\-ly tr^'-H-'' (u + v + l)-^du dv. 

Y Y 

The coefficient of u in the exponential in the second integral being now 
greater than that of v, we carry out first the integration in regard to u. The 
first part of K& is, as above, equal to 27r^. Thus 

Kj = 27r2—f { e-'^^^-^'^^-^^ (u + iy (v-{-iytr''---'^v~'' {u-i-v-i-iy^dudv. 

Y Y 

In extending the 'z;-contour of this last integral, instead of the ^/.-contour 
as in Ka, we add to it 

I du 6-^(i+^)-^^(2^ + l)^ {v-\-iyu-''-h~''(u-{-v+l)-^dudv 



= 27n [ du 6-(^-^)"/^ = — 47r2 ; 

Jy 



7i, as above, standing for a contour encircling — (1 + ^), and not zero. As 
before, the extended integral is zero, and therefore the unextended integral 
equals 47r^. Thus Kj = 27r^-— 47r^ = -— 27r^ for all values of h. 
We have 

Ltr = 00 I / (y) dy j j e-^y-^ a+i^) (^, ^_ l)r (^. _|_ ly ^^-r^-r (^^ ^^ + 1)-1 6^^6 ^v 

Jo •'v'y 

= -277/ (a;-) 
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and 

M' = 00 [ f{y) dy [ [ e-w2/-^(i+^) {u-{-iy {v-\-iy ir'^V'' (^6 + ^+ 1)~^ clu civ 

= 27ry(^ + ); 
and, therefore, finally 

n 

Ltn = a^tarCOrX = H/ (^ — ) +/(^ + )}> 


and if the function / is continuous, 

n 


It will be noticed that the success of the above method depended on the 
subject of integration in the beginning being a geometrical progression. 
If the general expansion for the unsymmetrical case (33) is treated similarly, 
the subject of integration is a double geometrical progression and may be 
treated similarly, but the analysis will not be carried out here. 

VII. — The Application of the (o-f unctions in the Theory of Frdbability. 

1. As was mentioned at the outset, the linear ^-functions have been 
recognised for some time as suitable functions for the purpose of representing 
a given frequency distribution to any desired degree of approximation. 
The method commonly adopted in fitting the coefficients is to make the 
successive moments of the required series equal to those of the given 
distribution.^' But it does not appear to have been noticed anywhere that 
this method is exactly equivalent to Hermite's method of obtaining the 
coefficients. This is an immediate consequence of the fact that the function 
e^'^hin is a polynomial. Thus the practical method by moments, which 
commonly appears as an artificial way of obtaining a series to represent a 
distribution of a particular type, is really the exact and unique way of 
representing an arbitrary function by a series of cw-functions. Edgeworthf 
has another method of obtaining a series of successive approximations to a 

r = t 

given distribution in the form exp {tkr { — cl j dxY j r !} {(2'7Th)~''e-^^''l^^'), which. 



r — s 



if the formal operator be developed, leads to a series of oo-functions. The 
method of obtaining the coefficients kr is that of identifying the moments, so 
that ultimately the series as obtained by Edge worth must agree with the 
direct expansion in w-functions. It might be, however, that, for statistical 
distributions, the method of development used by Edgeworth would give a 

"^ Charlier, loo. cit, 

+ ' Camb. Phil. Trans.,' Aug., 1904. 
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more rapid approximation. The argument as to the magnitudes of the 
coefficients kry however, appears equally applicable to the present series.: 
It is shown in the memoir quoted that the coefficients k^ rapidly diminish 
only if the standard deviation be small. Let this be assumed in forming 
the series anUnix), where it^ — (dldxy {{2ira^)~h~'^'^l^'''^], Then, if fin be the- 

r 

n\h moment of the distribution, and tt^ = I e ^'^IhO'du, we have 

J —00 



•00 

n 



fln^t{-y CtrG'^-'-lTn-rn 1 /(7^ - t)\ (LV)^, 


and in particular 1x2 = (27r)~^ {^oO"^7r2 + 2<227ro}, 

and hence, since a is the standard deviation, tt^ is zero. KTow, following 
Edge worth's argument, yu„ is of the order of magnitude cr^ and hence, by 
taking the above equations for different values of n, we see that the 
coefficient an is of the order cr^ — the word order being used in the same 
sense as in the work cited. Thus the coefficients tin will diminish with about 
the same rapidity as the quantities kn, while the formation of the equations. 
is much simpler. 

Pass now to the case of two dimensions. Eeference has been made to- 
Hermite's note on the expansion of a function in a series of the form 

which has been shown to be equivalent to an expansion in co-functions. It- 
is a problem of some importance to obtain a series to represent a frequency 
distribution in two correlated characteristics, such as, for example, is con- 
sidered by Perozzo in his memoir on the a^jes of husband and wife at 
marriage.* The lines of equal frequency given by Perozzo are far from being 
concentric ellipses, as they would be if the frequency surface were a Gaussian 
normal surface z = e~i (^^'+ 2/12:2/ +&2/2)^ j^ ig therefore worth enquiring whether 
a few terms of the expansion in co-functions can be made to give an 
approximation that is tolerably near. The theoretical possibility of repre- 
senting the distribution by a convergent series is not here considered, though 
the conditions of validity of the expansion are almost certainly satisfied. 
Hermite gives a method of obtaining the coefficients. Using the notation 

(X) = e-i{ax^+2hxy + bf-)^ ^ ^^^ _ ^~i^(bx^-2hxy + aij^)^ 

(d/dxr {d/dyycl> = (/) . U«, ,, (3/3^)- {d/dyy yjr = yjrY^n, 

TJni^n and Ymn are both polynomials of degrees (m + n) in x and y. lif(x, y) 
be expanded in the form S%a,n, n^^^m, n, it is proved that 

^^lUy n ^^ ^m, h\ \J \p^y) * m, n^X ay, 

•» — 00 J — 00 

^ ' Annali di Statistical 1883. 
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where a,^, ^ is a nnmerical constant depending on m and n alone. Now it 
will be mucli more convenient in practice to obtain the coefficients by means 



of equating the successive moments \ f (x, y) x'^^^dx dy to the corre- 

J — QO J —00 

spending moments of the series. This is equivalent to Hermite's process, 
V^», n being a polynomial ; so that a^, ^ above is a linear function of the 
moments. 

I Ur^'^yHx dy can be reduced by successive 

— 00 J —03 

integrations by parts. If m< r, we have 

11 Ur^^^^y^dx dy — {-'Y'ml y'^Ur-m, s dxdy = 0, 

whatever the values of s and ?^, and, similarly, the integral is zero it n<8, 
whatever the values of r and m. But if m = r, and n = s, 

iumn^^^'y^dx dy = (—Y'''^^ m I n I {\<f> dxdy, 
.and if m>r and 7i>s, 

1 [ur^'^yMx dy = (— y+^ (m {/(m—r) !) (n 1/(^—5) !) f^:^^~y "*& dy, 
und this last integral is again zero if m—r + ^i— 5 is an odd integer. 

Using the notation (l>rs = ipx'^y^dx dy and firs = fx^y^dx dy, if we 

multiply the equation /(^, y) = ^ars, ^hs by x^y^ and integrate, we have 

+ {r (r — 1) a2o(l>r-2, s + ^' - San^r-l, s-l + ^ (s— 1) %2#r. s-2) . . . , 

there being (r + 1) (s + 1) terms on the right, of which those are zero in 
which the sum of the subscripts of (j> is odd. 

Taking the first few equations, by choosing the origin at the mean of the 
distribution, and choosing the constants <x, &, h so that 

the constants aoi, aio, ago, «ii, ^02 all become zero. By this means the normal 
surface that best fits the distribution is obtained, viz., //'Oo<^/<^oo* 

But the fact that the contour lines in the example quoted are by no means 
concentric ellipses with centre at the mean shows that the fit is not yet 
.sufficiently good. 

The next approximation to the distribution is 

/Jm<p/4^00 + {/^BOUSO 4- 3/i2l'?^21 + 3/Xi2'?^12 + /i03^^03)/6^00- 
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The equations giving the coefficients of the fourth order terms are : 

Thus in terms of the moments firs the coefficients are quickly obtained^ 
and Edgeworth's argument as to the rapidity of approximation applies, 
equally well here if the standard deviations are sufficiently small. But the 
labour of obtaining the moments is very great, and a suitable example for 
testing the closeness of approximation is not available. What will be done 
here, therefore, is to show more or less generally how, in a simple case, contour 
lines can be obtained which are approximately of the type found in the 
example given above (see diagrams, Perozzo, loc, cit), 

2. An Example of Gontoitr Lines in Tivo Dimensions. — The example chosen 
here is that of ^ = g-M^^^ + y^)' and it will be assumed that there is symmetry 
about the line x = y. 

The most general first correction is a('i% + ^to3)4-/3('if2i + Wi2), where 
^30 = ^ (3ir — o(?\ 1121 = (t>y(l — ^% %2 = <^^ (1 — ?/2), toos = ^ (Sy — y^). 

In particular, if /3 = 0, and u is negative, ci(tcso + Uos) vanishes on the- 
straight line a? + y = 0, and the ellipse x^--'Xy-\'y^ = 3; on the positive 
side of x-\'y = it is positive outside the ellipse, and negative inside, and 
vice versd on the negative side oi X'\-y = 0. 

Thus the effect of the correction a (ttso + nos), while not altering the mean 
or second moments, is to introduce skewness about x + y = 0^ the maximum 
ordinate being displaced in the negative direction along x = y. Outside the 
ellipse x^'—wy'\'y^ = S the surface is depressed on the negative side of 
x-i-y = 0, and raised on the positive side. The zero-line of <^ + a (u^o + uos) is. 
shown in the figure (curve A) for a = — 1/10. 

The contour lines on the positive part of the surface will now be oval curves, 
surrounding the point of maximum ordinate with their greatest width parallel 
to X'i-y,diiid. crowding closely together in the neighbourhood of the zero- 
curve. These have not been shown in the figure to avoid complication. A 
variation in their form could be produced by introducing a terrii /3(%2 + ^2])?- 
but the effect of skewness is again the most noticeable. 

Consider now a fourth order correction, which for the end in view has been 
conveniently taken to be 

y(u22—iisi — ui3) = 7 {l'-a^'-y^ + 6xy—xy{x^—xy-^y^)}. 

The curve for which this is zero is drawn in the figure (curve B), and the 
portions of the plane are indicated in which it is positive and negative 
respectively. Thus the effect of this correction on the contour lines in the 
neighbourhood of the maximum ordinate is to contract them in the direction 



sso 



I he 0}-JbvMcU07is. 



/s = —y, and to elongate them in the direction x = y. A series of contour 
lines is drawn for 7 = 1/20; so that the lines shown are curves of equal 
Irequency for the function 

The curves are distinctly similar in type to Perozzo's curves, though no 
.attempt has been made at determining exactly what the coefficients should 
be. Probably a better approximation would be found with a rather larger 
coefficient for the third order correction, as the skewness is not so marked as 
it should be. 
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Figure diowing Contotir Lines of the Swface 

A— the curve ] 1 - ^V (8'/2^^ + 8%')} ^"^ ^""""^'^ = ^' 
the curve % = {d^/dx%^ - d^/dx^ - ^'^/^aB|/^ e-^ i^-^'J'^ = 0. 
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Similar curves can clearly be obtained where there is correlation to be 
taken into account. Tor by making the change of variable 

is turned into a normal surface with correlation, while c/dx and d/Bi/ are 
linear functions of d/d^ and 3/9/7, 
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(From tlie National Physical Laboratory.) 

(Abstract.) 

Attempts to prepare pure eutectic alloys of known constitution led to the 
discovery of discrepancies between the authors' experiments and the data 
on lead-tin alloys published by Eoberts-Austen. A complete redetermination 
of the equilibria of the lead-tin system was therefore undertaken, by both 
pyrometric and microscopical methods. Cooling-curves of the alloys taken by 
both inverse-rate and differential methods are given ; these, together with the 
microscopic data, lead to the equilibrium diagram shown in the figure. This 
differs principally from that given by Eoberts-Austen,* in that the eutectic 
point is placed at a concentration of 63 per cent, of tin instead of 69, that the 
eutectic line towards the lead end of the series terminates at a concentration 
close to 16 per cent, of tin, and that a series of transformations (along the 
line EFG) have been found in the solid alloys near the lead end of the series. 

The discrepancies of these results, as regards the eutectic composition, arise 
from the more delicate method employed in the present research. While alloys 
within 1 per cent, on either side of the true eutectic composition show no 
detectable difference in freezing- or melting-point, the presence of small 
amounts of excess of either constituent can be detected by the microscope, 
^nd in this manner the composition of the pure eutectic has been ascertained. 
As regards the solubility of tin in solid lead, it was found that the occurrence 
of the eutectic arrest-point in alloys containing less than 16 per cent, of tin 

■^ Eoberts-Austen, ^ Fifth Eeport to the Alloys Eesearch Committee of the Institutio 
of Mechanical Engineers,' 1897. 



